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In an attempt to simplify the design of separation filters in fre- 
quency-division multiplexed digital data transmission by fsk, it was 
found that the energy on one side of the fsk carrier can be reduced 
by suitably modulating the amplitude of the fsk signal. We present 
a simple technique for controlling the energy on one side of the 
carrier of an fsk signal. The method entails modulating the ampli- 
tude of the fm wave in unison with the frequency. In other words, if 
the energy in the region above the carrier must be reduced, the 
envelope is somewhat suppressed while the higher frequency is trans- 
mitted. The proposed hybrid modulation technique is applicable in 
frequency-division multiplexed digital data transmission by fsk and 
can often ease the design of separation filters. Here we investigate 
this modulation technique theoretically as well as experimentally. 
Our theoretical results include closed-form formulas for the various 
power spectra of interest. We also determine the optimum amount of 
amplitude modulation to achieve the greatest amount of sidelobe 
energy suppression for a fixed amount of total average power. The 
experimental investigations substantiate the analysis, and we make 
comparisons of the technique with alternative approaches. 

I. INTRODUCTION 

It has been known that complete elimination of one of the sidebands 
in fm can be achieved by modulating the envelope by a signal related 
to the baseband modulation. This scheme is known as ssb-fm. 1,2 In 
general, beneficial spectral shaping in fm is achieved only if the 
amplitude modulating signal is related in some manner to the fre- 
quency modulation. Also the amplitude modulation must be positive 
so that the zero-crossings of the original fm wave will not be altered. 

For random data, the spectrum of fsk is known to be symmetrical 
about the carrier frequency, and it appears intuitively clear that 
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reducing the carrier amplitude during the period while the higher 
frequency is being transmitted reduces the energy of the spectrum in 
this frequency range. 

We have observed experimentally that the amount of amplitude 
modulation critically affects the energy distribution. These preliminary 
experimental observations prompted us to examine this subject ana- 
lytically and to verify the analysis by further experimentation. 

Section II outlines the derivations of the am-fm spectral density 
when the baseband modulation is a synchronous digital data signal. 
Section III repeats the derivation for asynchronous data. Section IV 
discusses numerical and experimental results 

II. THE SPECTRUM OF AM-SYNCHRONOUS FSK 

To illustrate the principal ideas and to bring out the possible benefits 
of am-fsk, we use two simple models of digital fm and restrict ourselves 
to simple amplitude modulation mechanisms. We first outline the 
derivation of the spectrum for a general parameter set and then 
specialize to binary fm. 

Consider a stream of data symbols, oq, a u • • • , where each a n can 
assume real values independently and with equal probability. Now 
form the pulse train 

D(t) = £ a n g(t - nT), (1) 

where g(t) is a unit-height, T-second duration rectangular pulse. This 
signal will be used to modulate the frequency of a carrier wave. For 
the amplitude modulation, form another pulse train 

ao 

D(t)= £ b(a n )g(t-nT), (2) 

where b(x) is a nonnegative function. Next, apply (1) to an fm modu- 
lator and multiply the output by (2). These operations result in an am- 
fm wave which is represented by 



S(t) = 23(*)expj i 



I D(t')dt' + \, 



u c t + ud D(t') dt' + d 



(3) 



where w f is the carrier angular frequency, co</ is a constant of propor- 
tionality, and 6 is a uniformly distributed phase angle on the interval 
[— 7r, ir\ We use the complex representation of fm for calculational 
economy without loss of generality. 

Because D(t) is always positive, D(t) can be recovered by conven- 
tional fm demodulation methods, 3 and the role of the amplitude 
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modulation here is only to shape the spectrum of the fm wave. The 
manner in which spectral shaping is accomplished is our main interest 
and motivates our investigation of the spectral properties of (3). 

Denote the spectral density of S(t) by G(co). A convenient definition 
of the spectrum is the limit 



\jm±E{\Z(o>,N)\ 2 }, 



(4) 



where 



Z(a, N) = 



S(t)e-^ 1 dt 



and where E{-} denotes the ensemble average. 

The sequence of squared magnitudes of Z(u, N) is calculated in a 
straightforward manner: 



\Z{u,N)\ 2 = 



J-{k+\)T r 

expj i 
kT I 



N-\ ■■■•'' Il7 

£ b{a k ) | ex pi 

k-o 



(co c - io)t + to d £ a n 

n-0 
t-nT 



By recognizing that the phase angle, 

\p(t) = (w - u c )t - io<t X a n 



g(t') dt' 



g(t') dt', 



(5) 



must be a linear function of time in the interval kT < t < (k + \)T, it 
is possible to put (5) into the form 



1 2 _ v t/_ \t/_ \jMs-k): 



\Z(u,N)\ 2 = I b(ak)b(a.)tf" l - k) F(p-a k u d )F*(p-a.a d ) 

k,s-0 

{r*-i s-i 

ia J] a„ - J] a„ 

where a = co^T, i» — w — to c , * is the complex conjugate, and 



(6) 



F(x) = 



,-«f 



* 



[See Ref. 4 for detailed manipulations leading to (6).] 

To get the spectral density, (6) is first averaged with respect to the 
identically distributed random variables, ao, fli, «2> • • • , and then the 
limit in (4) evaluated. Manipulations similar to those carried out in 
Ref. 4 lead to the desired result, namely, 
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G(to) = lim-£|Z(<o, N)\'- 



= E a 



b(a)Si 



vT - aa 



+ 2 Re 



_-i>r 



1 - e- ivT C(a) ] 



':"' , (r'l ota\_. . 
\b(a) 



e^Si 






(7) 



where a is a data symbol, C(a) = EJiff—), and S,-(x) = sin x/x. The 
above result for G(w) holds only when | C(a) \ < 1. When | C(a) \ = 1, 
the spectrum contains delta functions. These can be determined with- 
out too much difficulty, 4 but, since in this application | C \ < 1, we do 
not need the most general results. 

In our application, the data are binary so that a = ±1. We also 
specialize to a = u d T = tt/2. For this special case, the spectral density 
formula reduces to a very simple form. 

Let 

6(1)= A 

and 

b(-l) = B. 

Now evaluate explicitly the averages specified in (7). 



+ Re\e 



-i7V 






Making use of the identities 



2 ,'fL^W^I-l-s»,T, 



2l'^IS? 



2 
a+Tv 



eq. (8) becomes 



G(p) = K 



= 1 + sin pT, 
(a = tt/2), 

B 



cos 2 i>T, 



a-Tv a + Tv 
where K is an unimportant normalization constant. 
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(9) 



Several checks reveal the physical reasonableness of (9). When 
A = B, we have no amplitude modulation and hence S(t) becomes a 
constant-envelope, continuous-phase fm wave. In this case, 

r i 

G(p) 



JvTY-c? 

which is a well-known result. 5 

Notice that, in this case, because of phase continuity, the spectrum 
behaves like 1/to 4 for large to. When A # B, the spectrum behaves only 
as 1/co 2 , as o> becomes large. In the extreme situation when A = 0, we 
have 100-percent amplitude modulation and, in fact, no frequency 
modulation. The only tone transmitted is at u c -u d , and S(t) is an on- 
offsine wave with frequency w t - u d . The spectrum of this wave is the 
well-known (sin x/x) 2 . Put A = in (9) and, with a = u d T = tt/2, obtain 

Tsintto - (w c - u d ))T'\' ii 

G(io) <*> 

(u - (u c - u d ))T 

which is in agreement with physical intuition. 

III. THE SPECTRUM OF AM-ASYNCHRONOUS FSK 

The baseband asynchronous fsk wave is constructed as follows. On 
the interval [0, T], pick a set of points at random and arrange them 
such that 

= to < k < fe < • • • t N m T. (10) 

Define a set of functions 

g *" (t U "{0, otherwise, (11) 

where 

A„ = tn+i ~ t n . (12) 

In terms of (10) and (11), construct the baseband signal X(t) as the 
following time series 

N-l 

X(t) = 2 a n g An (t-t n ), (13) 



where a =' (ao, flti, •••, a^-i) is an arbitrary set of identically 
distributed random variables. 

The instantaneous phase \p(t) in this case is 

ifr(t) = uct + aid I x(t') dt' + 6 
Jo 
= «M*) + 0, 0<^<T. (14) 
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The constants w c and u d have the same definition as in the previous 

section. 

As before, let b{x) be a nonnegative function and construct the am- 
asynchronous fm wave as follows: 



S(t)= £ b(an)gA m (t - t n )e 



i+i(t)+i6 



(15) 



Now define the truncated Fourier transform 

Z(o>, T) = I S(t)e~ iul dt 
Jo 

= e w Yb(a k ) J apptyift) - «0] A 



A=0 



"- 1 b(a k ) 



= e" % -4 



'o IAa 



where 



exp{ i £ A„A„ [ - expj i £ A„A„ 



An = COrfOn — W + C0 C . 



(16) 



To proceed further, the magnitude squared of (16) must be averaged 
with respect to the set of random vectors 

A = (Ai, A 2 , •••, An) and A = (Ai, A 2 , • • • , \n). 

In (15), the A„'s are the intervals between transitions. We shall assume 
that these intervals are independent and identically distributed. On 
the other hand, the random variables A„, which are linear functions of 
the amplitudes a n , are not independent if only actual transitions are 
considered. Clearly, for observable transitions a n i* a n -i and therefore 
adjacent amplitudes cannot be independent. 

This difficulty in the analysis has been overcome in Ref. 6 where 
asynchronous fm without amplitude modulation is treated. Our prob- 
lem here is a slight modification of the pure fm case, and identical 
reasoning and techniques to those exhibited in that reference can be 
used to obtain a closed-form formula for the power spectrum with 
amplitude modulation. Since the algebra leading from (16) to the final 
result is rather involved, we shall not repeat it here and will give only 
the final formula. 

Thus, for the special case where the number of t points in a fixed 
interval T obeys the Poisson probability law with parameter v = 
(l/TV/i = average number of real transition per second, we get for the 
power spectrum 
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where 



and 



G(u) = lim-£|z(to, T)| 2 

T-.x. 1 

_ [AT(J3 + u d ) - BT(u - u> d )f 
(n>) 2 + (T 2 & 2 - <o 2 d )) 2 



T = \/2v, to = w — w c 



6(1)= A, b(-l) = B. 



(17) 



Our experimental work is based on synchronous data, and formula 
(17) is included only for completeness. Clearly, the spectrum for the 
asynchronous data exhibits similar properties to that of the synchro- 
nous case. 

IV. NUMERICAL AND EXPERIMENTAL RESULTS 

In Figs. 1 and 2, we have plotted several curves exhibiting the power 
spectrum G(f) [eq. (9)] for different values of the parameter 




Fig. 1 — Power spectral densities for fsk and am-fsk. 
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Fig. 2 — Power spectral densities for am-fsk. 

p = A/B.f An fsk transmitter with amplitude modulation (am-fsk) 
has also been implemented. The measured power spectra for different 
values of p are shown in Figs. 3 and 4. In the experiments, we have 
chosen a carrier frequency of 250 Hz and a bit rate of 75 b/s. The 
agreement between the analytical and experimental results is seen to 
be excellent. For p = 1, we have a conventional binary fsk signal. From 
Figs. 3 and 4, we see that when p decreases the amount of power in the 
right sidelobes decreases at first and then increases again. At the same 
time, the amount of power in the left sidelobes increases steadily but 
at a slower rate. Quantitative results are given in Fig. 5, where we have 
plotted the residual powers P u - S/>t G(f) df and Pl = -\-?,t G(f) df 
in the right and left sidelobes, respectively, as a function of the 
parameter p. For the two sets of curves, corresponding to two different 
values of f/T, the minimum value of the residual power in the right 
lobes is obtained for a value of p approximately equal to 0.75. It is 
possible to derive a closed-form expression for the value of p which 
gives the minimum residual power as a function of far. The derivation 
is given in the appendix. For fa = 1, the optimum theoretical value of 
p is 0.7623. 

am-fsk has potential applications in situations where the power in 
the right or left lobes of a conventional fsk signal might be detrimental 



f In the following we always assume p < 1. The spectra corresponding to p > 1 are 
simply the mirror images with respect to the carrier frequency of the spectra correspond- 
ing to p < 1. 
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to the proper functioning of a contiguous neighbor. Filter requirements 
can also be eased if power in the interfering lobes can be reduced. 
Clearly, methods other than amplitude modulation exist to reduce the 
power in the lobes of an fsk signal. One obvious method is simply to 
filter some of the power. However, in a digital environment, such an 
approach can be rather costly. Typically, a good nonrecursive filter 
requires about two orders of magnitude more computation time than 
is required for the implementation of the amplitude modulation de- 
scribed here. A second possibility to decrease the amount of interfer- 
ence produced by the lobes of an fsk signal on another signal consists 
in moving away the carrier frequency of the fsk signal. To investigate 
this possiblity, we have computed the amount A f c - T by which the fsk 
spectrum should be shifted with respect to the am-fsk spectrum in 




p = 1 




P = 0.75 



Fig. 3 — Measured power spectra of fsk and am-fsk (10-dB per division along the 
vertical axis). 
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p = 0.5 




p = 0.25 



p = 
Fig. 4 — Measured power spectra of am-fsk signals. 

order to get the same residual power above the point f/T = l.f The 
results are shown in Fig. 6. In the neighborhood of the point p = 0.75, 
the shifting A/ c of the carrier frequency has to be close to 1/T to get 
the same performance for the fsk and am-fsk signals. 

Finally, a third way to reduce the power in the sidelobes of the fsk 
signal consists in simply decreasing the total power of the signal. In 
Fig. 6, the dashed curve shows by what amount the power of the fsk 
signal must be reduced with respect to the power of the am-fsk signal 
in order to have the same residual power beyond the point fJT = 1. It 
can be seen that around the optimum point, near p = 0J5, this is 
achieved when the power of the fsk signal is about one-tenth the 
power of the am-fsk signal. Figure 7 gives a more detailed comparison 



f The point f/T = 1 corresponds to the center of the first lobe of the nonshifted signal. 
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between the residual power for fsk and am-fsk as a function of f/T in 
the case where p = 0.75. 



APPENDIX 

The Optimum Value of p 

The residual power in the upper sidelobes of the power spectrum is 



given by: 



P„ = 



A B 

+ 



cl-Tv a+ Tv 



cos 2 Tv dv, 



(18) 



where a = •n/2 and v f > ot/T. The functional P u is minimized under the 
power constraint 



P = A 2 + B 2 = Constant. 



(19) 



This variational problem can be solved in a straightforward manner 
by introducing a Lagrange multiplier and evaluating 



40 




Fig. 5 — Residual power in the upper and lower side lobes of am-fsk. 
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Fig. 6 — Curves showing by how much the fsk signal should be shifted or attenuated 
in order to have the same residual power as am-fsk. 



dA dA 
dB dB 



(20) 
(21) 



Replacing P„ and P by their expressions and letting p = A/B, we have 



a-Tv 

1 

a + Tv 



" +. ' 



a- Tv a + Tv 

o 1 

+ 



- Tv a+Tv 



cos 2 pT dv + Xp = (22) 



coflfrT dr + A = 0, (23) 



Eliminating X between (22) and (23) and redefining some variables, we 
get 



Kz — Kz 



P ~P 



K, 



-1 = 0, 



(24) 
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where 






dx 



K , = I J^dr 



y,T 



(a - xY 



COS X 

Kz = I — dx. 



*,t 



(a + xY 
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Fig. 7 — Residual power in the upper sidelobes of fsk and am-fsk. 
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These integrals can be expressed in terms of tabulated sine and cosine 
integrals. 

After some algebra, we get 

Kl = I in 2veT ~ m + - ci{2v,T - 77) - - dfr/T + 77) 

77 2v/T + 77 77 77 

K 2 = — J (1 + cos 2p,T) - 2 s\(2v/T - 77) 

2v/r — 77 

K 3 = — J (1 + cos 2v/T) - 2 si(2v/T + 77), 

2v/T + 77 



where 



and 



, sin x 
si(2) = - ax 



1 cos x j 
ci(z) = - ax. 



Tables for the si(z) and ci(z) functions can be found in Ref. 7. 

The optimal values of p are the roots of the quadratic equation (24). 
The root that is positive and less than one is the solution to our 
problem. It is given by 



Popt — 



K.2 — K$ If K.i — K3 



2Kx ' M\ 2Ki 



+ \l ' ) +1. (25) 
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